function outcor = f_heatmap_subplot(Nr,Nc)
a = 0.01;
b = 0.11;
cc = 0;
Lx = (1-a*(Nc+1))./Nc;
Ly = (1-b*(Nr+1))./Nr;
for i=Nr:(-1):1
    for j = 1:Nc
        cc = cc+1;
        outcor(cc,:) = [a+(j-1)*(Lx+a) b+(i-1)*(Ly+b) Lx Ly];
    end
end
return